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We discuss a charge qubit consisting of a circular array of Josephson junctions. The two-level system 
we consider couples the two charge states through a higher order tunneling process thus making it 
possible to achieve a long relaxation time. Using the spin-boson Hamiltonian, we estimate decoher- 
ence due to ohmic as well as 1// noise. We simulate the quantum mechanical measurement process 
by studing the density matrix of the qubit and a capacitively coupled single-electron transistor that 
measures the charge. 
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I. INTRODUCTION 

In recent years much effort has been spent on the search for quantum two-level systems, qubits, that can be coherently 
'■^ ■ controlled long enough for a sequence of controlled unitary operations to be performed on them. The ultimate goal is 
to build a quantum computer out of these qubits. Proposals for qubits based on a variety of physical systems exist, 
each with its pros and cons. Here we consider solid state charge qubits based on Josephson junctions (JJ) arrays. 
These have the advantage of being relatively easy to manipulate and the prospects for large scale manufacturing are 
\ comparatively good. However, they suffer from severe decoherence effects. The single Cooper-pair box (SCB) is the 
C]3 ■ simplest proposal for a charge based JJ qubiti For the SCB, superposition of charge states was observed by Bouchiat 
g \ et al£ and coherent evolution was demonstrated by Nakamura et alS^ 

A generalisation of the SCB, which we call the circular array (CA), was introduced by Schollmann et alX This 
circuit consists of an array of tunable JJs in a circular geometry. The CA is similar to the SCB and many results can 
be taken over mutatis mutandis. The main difference is that the two charge states of the qubit are coupled through 
£j ' a higher order tunneling process. Turning the coupling off then allows the tunneling rate to be made very small - 
leading to slow relaxation, and a long time to perform the measurement. This is the key element of the quantum 
sample and hold (QUASH) measurement strategy^ 

In this article we perform a more detailed study of the circular array. In particular, we consider the effect of voltage 
fluctuations in the circuit (ohmic noise) as well as 1// noise, believed to be caused by background charge fluctuations, 
j— ( ' and calculate the relaxation and dephasing times for these types of noise - extending the previous treatment^ We 
. also study the measurement of the qubit's charge by a single-electron transistor (SET) coupled capacitively to the 
CA. This is done by numerically determining the time development of the density matrix following the treatment of 
O ; Makhlin et al. 6 for the SCB. 
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II. THE CIRCULAR ARRAY 



I ■ The circular array consists of two arrays with N identical JJs each - these arrays are connected in series and separated 
by a capacitor Co, thus forming a circular geometry. Each JJ, which is a small SQUID, has capacitance Cj and a 
Josephson energy Ej — Ej($(t)) which can be tuned by altering the magnetic flux $ through the SQUID loop. 
The lead connecting the two arrays is grounded to allow charge to tunnel in and out of the circuit. There are 2N 
small islands, i — 1,2,.. . 2N, each characterized by the number of excess Cooper pairs rij and the phase of the 
superconducting order parameter (f>f, these are quantum mechanically conjugate variables: [<f>i, nj] — iSij. Each island 
charge is externally controlled by a gate voltage Vi(t), applied via a small capacitor C g . Fig. ^ shows the circular 
$Jj ■ array together with the SET that measures the charge on one of the islands next to Co- 

A qubit should have two states separated by a large gap, SE, from higher energy states and be weekly coupled to 
the environment to avoid rapid decoherence. The CA fulfills this if Co ~ Cj C g . The energy scales present in the 
system is the charging energy for a Cooper pair Ec = (2e) 2 /2Cj ~ SE, the Josephson energy Ej, the superconducting 
gap A and the temperature fc^T. In order to avoid quasiparticles in the system at low temperatures, the qubit is 
constructed so that A is the largest energy in the problem. Furthermore, we choose the qubit to be in the charge 
regime Ec 3> Ej, and impose Ec 3> &bT to avoid thermal excitation of higher energy charge states, thus: 



A> £ c > Ej,k B T 



(1) 



2 




FIG. 1: The circular array with a SET connected to island N. The box symbols (without cross) stand for normal junctions. 
Island N on the CA and the island of the SET is marked with dotted boxes. 



The Hamiltonian of the CA is 

2N 

H = H C +Hj = -J2 QiC^Qj - Ej cos(& - , (2) 

*»J=1 (*i) 

where Qi — 2e(rii — n g ^) is the effective charge on island i - here n g ^ = C g Vi/2e is the gate charge on the island. 
Cij is the capacitance matrix - its nonzero elements are: Cn,n = Cn+i,n+i = Co + Cj + C g , Ci.i = 2Cj + C g , 
Cv+i.n = Cn,n+i = —Co and C + i^ = C.i +1 = — Cj, where i ^ N,N + 1. The matrix is symmetric and 

= C2N +i-i,2N the inverse matrix C _1 has the same symmetries. The sum over the Josephson terms in 
is taken over all pairs of islands connected by tunnel junctions. 

Since we are studying a charge qubit it is convenient to write the Hamiltonian in the charge basis |n) = 
\niri2...n2N)- The charging energy term simply becomes He — EcCjJ2n( n ~ n £/)*C _1 ( n — n s )|n)(n| where 
rig = (n g i,n g 2,---,n g ,2N)i and using that \m) = J** tyie-^'lfa), which holds since n, and fa are conjugate, 
the Josephson term becomes 

Hj = -^-J2 II + l )\ n i - !) + l n « - l )\ n i + !))( n l ( 3 ) 

and the total Hamiltonian is H = Hc(V) + Hj($>), where we have indicated the dependence on the external control 
parameters V = (Vi, V2, ■ ■ ■ , V-jjv) and $. 

The two-level system - ie the qubit - we consider consists of states |J.) with one excess Cooper pair on either of 
the islands N, N+ 1 neighbouring Co: |T) = |0 ... 10 ... 0) and |J.) = |0 ... 01 ... 0), where the ones are for island N and 
N + 1 respectively. For Ej — these two states are degenerate if n g jy — n g jy + i — | and n g i = for i 7^ N, N + 1. 
If, in addition, Co ~ Cj ^S> C g , then the energy gap to higher charge states is SE ~ Eq. Restricting ourselves to 
a finite charge space, the Hamiltonian H can be diagonalised numerically. In Fig. [21 we show the energy spectrum 
for the N — 2 CA as a function of n 3i 2- The other parameters are n g ^ = 1/2, Co = Cj = 100C g , Ej = 0.2Ec and 
—2 < n% < 2 (this restriction gives a negligible error). For n g> % ~ 1/2, the two lowest energy levels (which are linear 
combinations of the states ||), \\) with an excess Cooper pair on island N or N + 1 respectively) form a two-level 
system with a large gap to the higher energy states. 

We conclude that when Ec S> Ej, Cq ~ Cj 3> C g and n ff w (0, . . . , |, |, . . . , 0), it is a good approximation 
to restrict the Hubert space of the Hamiltonian in Eq. (J2J to the states ||) and ||) defined above. We write the 
Hamiltonian of this two- level system in spin- 1/2 notation 

H ct ri(t) = ~B z {5V{t))a z - ^B x (<5>(t))a x , (4) 

where er^ are the Pauli matrices in the basis |J,)). This Hamiltonian controls the qubit - unitary operations 
can be performed on the qubit by tuning B z and B x via the external parameters SV = V/v+i — Vpj and $. Writing 
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= (0, 



l + Srig l — 5tlg 

2 ' 2 

n- 1 



t) = 2AjqEc&n g , where 



,0) where Sn g = C g SV/2e, we find B z = (|| H c II) - <T \H C 
An = Cj(C NN — C N 1 N+1 ). An can be calculated numerically for given capacitance matrix, however, we can also 
perform an expansion in C g /Co which is valid as long as N is not too large. Using Cramer's rule for the elements in 
C , we express An in terms of cofactors and expand. Assuming, for simplicity, Cj = Co, this gives 



B, 



2N 



2N+1 



C g N + l 
C 6 



+ O((C g /C ) 2 ) 



x E c Sn g 



(5) 



(For N = 2 and to leading order in C g , this reproduces a previous result^) B x gives the cotunneling rate of a Cooper 
pair from island N to N + 1 via the 2N junctions and and its leading contribution is obtained by (2N) th order 
perturbation theory, hence B x ~ (Ej / Ec) 2N ■ The exact numerical factor is not very illuminating 

determined for not too large N. For the N = 2 CA, with parameters as in Fig. [21 we have B z T 



it can however be 



1 _ l£<L 

2 Co 



E c Sn n = 



0.796E c Sn g and B x = ^{Ej/E c yE c = 0.0167^c*- (A numerical diagonalisation gives B z = 0.7960E c Sn 9 and 
B x = 0.0132£'c' _ in reasonable agreement with the expansions taking into account that the ignored terms are of 
order O{{C g /C ) 2 ) ~ 10~- 4 and 0{{Ej/E c f) - 1(T 4 respectively.) 



Diagonalizing (J2J gives the eigenvalues ±^AE where AE = \JB%. + B%, with corresponding eigenvectors 



cos : 



T) 



I), and 
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I), where 9 = aicta,n(B x / B z ). If r denotes the Pauli matrices in 



the energy eigenbasis, then H ctr i can be written in the compact form 

Hctrl = —-^AEt z . 



(6) 



III. ELECTROMAGNETIC NOISE 



The practical usefulness of a circuit like the CA as a qubit is ultimately limited by the coupling to external degrees of 
freedom. These lead to decoherence of the qubit state and hence loss of quantum information to the environment. The 
generic behaviour of the evolution of the qubit depends on the strength of the coupling and one identifies two regimes: 
The "Hamiltonian-dominated" , where the coupling to the environment is weak enough for the time evolution of the 
qubit to be governed by the qubit Hamiltonian H ctr i , and the "environment-dominated" , where the coupling to the 
environment is so strong that it determines the dynamics of the qubit. In this article we consider only the Hamiltonian- 
dominated regime. Note, however, that even if the coupling to the environment is weak under normal operation of 
the qubit, it becomes environment-dominated if H ctr i ~ 0, which may happen during the qubit manipulations. 

In the Hamiltonian-dominated regime the evolution of the qubit is conveniently described in the energy eigenbasis 
(|— ), |+)). The interaction with the environment leads to a decay of the off-diagonal elements in the qubit 's density 
matrix with a charecteristic time t v , the dephasing time, 



(r±(i)) = (r ± (0))e^ AB * e ^ ; 



(7) 
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whereas the diagonal elements of the density matrix decay to their thermal equilibrium values with a characteristic 



time T re i ax , the relaxation time, 



(r.(t)) = r z (oo) + (r 2 (0) - r z (oo))e^ 



(8) 



where the thermal equilibrium value is t z (oo) = tanh(A£ , /2fcsT). 

A Josephson junction charge qubit is sensitive to various electromagnetic fluctuations in the circuit; we follow 
standard practice and model these with the "spin-boson" modeli with an Ohmic spectrum. In addition to the noise 
caused by these fluctuations one observes 1// noise, which is believed to be due to background charge fluctuations in 
the substrate. Following Shnirman et al^ we model this phenomenologically using again the spin-boson model but 
now with a l//-spectrum. 

The spin-boson model describing the qubit interacting with the environment has one independent bath for each 
island in the CA: 
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SB 



H, 



ctrl 
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E^ 



(9) 
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7n a iUj a iX~ 



is a bath of harmonic oscillators with 



Here, H ctr i is the qubit Hamiltonian (@J, H l B - , 

coordinates x a %, momenta p a i, masses m a i and frequencies Lo a i- The baths lead to voltage fluctuations Xi = ^ a C a 
coupling to a z . (C a i is the strength of the coupling between the qubit and the a'th oscillator in bath i.) 

The effect of the environment is completely characterized by a spectral function, which for the spin-boson model 
in Eq. has the form 



j^) = |E— ^ 

2 ^— ' ni,,,~„, 



i) = \ E a ° h -' 
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LU S &(uJ c — U>) 



(10) 



To obtain the second equality one assumes that J(w) can be written as a power of lu up to some cut-off frequency 
lo c which is assumed to be large compared to all other frequencies in the problem. The parameter s reflects the 
qualitative nature of the environment, and a s = £\ a l s is a dimensionless measure of the strength of the coupling. To 
maintain a s dimensionless for all s, an additional frequency scale woo enters for s ^ 1. 

We model the voltage fluctuations SVi on island i by adding an impedance Zi(ui) in series with Vi, see, eg, Ingold 
and Nazarov^ or Makhlin et al£ This impedance then has Johnson-Nyquist fluctuations SVi between its terminals 

where Zn(u>) is the 



that in the spin-boson formalism correspond to the spectral function Ji(u>) = ujRe Zuito) 



where C Sj 



total impedance between the terminals of Z^. From Fig. [I]one finds Zn(u>) 

Cg 1 + (siCo)~ 1 , with Si = ?' _1 + (2iV+l — i)^ 1 (assuming Cj = Cq^S> C g ). Following standard practice, we assume that 
the noise is purely resistive, ZiiuS) = and if furthermore Ri <C l/uiC^i (which holds for realistic Ri and C g since 
Cs^ < C g ), we obtain Ji{uS) = toRi, which is linear in uj and hence corresponds to s = 1 in Eq. (|10|l . Using this we can 



obtain the total spectral function J(u>) for the circuit. From Eq. JjJJ we identify X 



■\8B Z 



where SB, 



(I \H C \ i 



> - (Tl H C |t) 



Simplifying this expression yields X — 



-eg^ Eifi AW, with At = Cj[C 



-l _ f-t-i 

N+l,i °JV,j 



The fluctuations SVi are expressed in the oscillator coordinates, SVi = Y]„ bi, 



form of Eq. JUJ, Jj(o 



-T ■ 

2 t—ia 



-5(U> - Ldia) = uRi 



This gives X = -e^E^i^Ea^ 



with the spectral function in the 
— - / ; | . •,„.'',„ and hence, 



J(w)=W E s : 



9 A 



Ri. Comparing this to Eq. I|10|l . with s = 1, we find 



(11) 



where Rk — hje 2 — 25.8 kfl is the quantum of resistance. Assuming the islands to be nearly identical, it is reasonable 
that Ri are approximately the same for all islands. Setting i?j = R, i = 1, . . . 2N, we are left with the factor E; Af, 
which we can calculate numerically for given capacitance matrix or as an expansion in C' g /Co. The expansion is 
analogous to the one for B z in Eq. (JSJ - assuming for simplicity Cj = Cq we obtain 



AN(N + 1) 
3(2N+l) 



2C g 
15C 



(N 2 + N + 3) + O((C g /C ) 2 ) 



Rk \ Co 



(12) 
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(For N — 2 and to leading order in C g , Eq. (|12|1 gives a\ — §^"(§jl as derived previous. 4 ) Fig. [31shows ai as a 

function of N for different C g /C - The numerical result is shown as crosses, and the expansion to order (C g /C ) 3 
as squares. For C 9 /Cq = 1CP 2 and 10~ 3 , the expansion agrees very well with the numerical result (at least when 
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FIG. 3: ai as a function of N for different C g /Co- The numerical result is shown as crosses and the analytic expansion as 
squares, ai is given in units of R/Rk and the scale on the y-axis is logarithmic; Ri — R for all islands, and Cj = Co- 

N < 10). For Cg/Co — 0.1 and 1, the expansion becomes negative at N > 8 and N > 2 respectively. Therefore no 
squares are seen in this region. In a typical circuit we may have R ~ 50 ft, yielding R/Rk ~ 10~ 3 . For C g <C Co, 
oil <Cl for all realistic N. 

The characteristic times T re i ax and r^, can be calculated within the spin-boson model using perturbation theory or 
path integral methods^ For s = 1 they are 

AE , AE . 9 „ 



relax 



^2 sin 2 ■ S x {AE/h) = vrai— - coth — — 



sin 



— 2 relax 



1 . , 1 2fc B T 2 

-2 • 5x(0) = -T re i ax + itoti6 — - — cos i 



where 



Sx(uj) = ([X(t),X(t')]+) u - 2?iJ(w)coth 



huj 

2k B T 



(13) 
(14) 

(15) 



is the Fourier transform of the symmetrized correlation function at thermal equilibrium. These results can be applied to 



-cos 2 9T 



the CA simply by substituting a.\ from Eq. 1)11(1 . One defines the pure dephasing rate T* , as T v — |r r 
The Hamiltonian-dominated regime is realized when AE 3> TiT* v (at least for s — 0, l 12 ). For ohmic damping, 
T* = 2naikBT/h, and the condition becomes AE 3> aiksT. 

Noise with a power spectrum proportional to the inverse of the frequency is observed in many physical systems. 
This 1/ / noise is seen also in J J circuits^ 1 -, where it is believed to be caused by background charge fluctuations^ 3 Here 
we follow Shnirman et al,® and model this noise with the spin-boson model. For s — and u> <C ksT/h, (|15fl gives 
the wanted power spectrum: Sx(to) = E^^/oj with — 2Trhctou>ookBT (here T is an adjustable parameter). Sub- 
ohmic environments (ie those for which < s < 1) have not been much studied as it was believed that they rapidly 
localized the system in one of the oveigenstates for any strength of the damping. However, it has been realized that 
this is true only for large damping, whereas for weak damping the system behaves coherentlyi*. The times T re i ax and 
t v have been calculated for 1/f noise for the single bath spin-boson model in the Hamiltonian-dominated regime^. 



The relaxation rate is 



hAE 



(16) 



The dephasing rates are only known for 6 — 0,7r/2; they are 



E 



i// 



7T TlUJir 



for 9 = 



(17) 
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and 

= St; tt- In t^tt; l'< »■ " = » ' -> • i Is I 

v hAE2n hu ir AE 1 K ' 

both with logarithmic accuracy in E\i^. Here, tu ir is an infrared cut-off frequency which can be experimentally 
determined. To determine the times we need to determine Ei/jw& 



IV. THE QUANTUM MEASUREMENT OF CHARGE 



The state of the qubit is inferred by performing a quantum measurement of the charge using a SET coupled capacitively 
to island N of the circular array see Fig.^ We simulate this quantum measurement by studying the time development 
of the density matrix describing the CA and the SET - we follow closely Makhlin et alfi where the corresponding 
problem is treated for the SCB. 

The SET is a circuit with one normal island surrounded by two junctions connected to normal electrodes and a 
capacitor Cg T , see Fig.^ During manipulations of the qubit, no current flows through the SET - the SET is turned 
off - this is achieved by setting the transport voltage to zero, Vtr = 0, and tuning the gate voltage Vg ET away from a 
degeneracy point so that the Coulomb blockade suppresses the tunneling through the SET. When the manipulations 
of the qubit arc done and one wants to read out the result, then the Coulomb blockade is turned off by tuning the gate 
voltage to a degeneracy point and a transport voltage is turned on leading to a tunneling current through the SET. 
This current depends on the state of the qubit through the charge on the capacitor Ci nt . This leads to a measurement 
of the charge of the qubit. 

The density matrix for the CA and the SET can be written as p = p(i, i', M, M', m, m')(t) after one has traced out 
the microscopic degrees of freedom in the left and right electrodes and in the island of the SET. Here, i labels the 
state of the qubit, M is the number of (excess) electrons on the island in the SET and m is the number of electrons 
that have tunneled through the SET. It is possible to derive a master equation for p as an expansion in the SET 
tunneling terms>i^ For low temperature and small Vt r , only transitions between two adjacent charge states of the 
SET need to be taken into account (these states are assumed to be M — 0, 1 below). If one furthermore assumes 
that the tunneling is instantaneous then one obtains a set of simple equations for the diagonal matrix elements 
p^(rn,i) = p(i,j;M,M;m,m)(t). In terms of the Fourier transformed quantity p l ^(k,t) — J2 m e ~ lkm PM( m , t), the 
final form is a system of eight coupled differential equations (for each k) 

n ±(h\,( %[HctriM \ = (- f t e- ik r R \(p \ 
dtyp-J \i[H cM + &H intl p 1 ]) \t L -t R J\pJ ' 1 > 

where pu is short for the 2x2 matrix p^(k,t). H ctr i is the qubit Hamiltonian with a renormalised capacitance 
matrix (Cn,n - * Cjv.at + Cmt) due to the presence of the SET, SHi nt — Ei nt a z is the coupling energy where Ei nt is 
determined by the capacitances. The tunneling rates T^/ R are 

FlPo = Fl/5 + ira L [5H int , p ] + , 

fflPa = r fl p! - na R [5H int , px] + , (20) 

where «l/_r and ^ l/r is the tunneling conductance and the tunneling rate for the left/right junction of the SET. 
We want to study the current through the SET as a function of time - this is obtained from 

P(m, t) = Y, P(*. h M, M, m, m)(f) , (21) 

i,M 

which can be interpreted as the probability that m electrons have tunneled through the SET during time t. To 
compute P(m,t), we solve the differential equation (|19|) with suitable initial conditions. We assume that the qubit 
and the SET are disentangled initially, p(t = 0) = p^ ® Pq ET , and that the qubit is prepared in some state \ip) — 
a IT) + P Pq = I^K^I- At time t = 0, no electrons have tunneled through the SET, thus P(m,0) = 6 m £ and 
P{k, 0) = 1 for the Fourier transform. From the definition l|21|) . we find P(k, 0) = Tr(po) +TV(pi) and we may choose 
(p ,Pi)|t=o = (0,p[f ). 

To calculate P(m,t) for the N = 2 CA we need values for the parameters E int , a^, a R: and T R for the SET, 
and B z and B x for the qubit. We choose C g = C mt = C^ ET = \C T = j^Cj = -^C = 0.032 fF, a L = a R = 0.03, 
Tl = 30 peV and T R — 100 peV (the qualitative result will not be too sensitive to this choice). Hence Eq — 100 peV, 
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FIG. 4: P{m,t) as a function of m for ten different times t. In (a), B x — 0, B z = 9.0 peV and \a\ 2 = 0.75, whereas in (b), 
B x = 1.3 /ueV, B z = 15 ^eV and \a\ 2 = 0.50. 



£i„t =6.1 /xeV and B z = 808n g /ieV. At the start of the measurement B x is turned off, B x ~ 0. B z is kept close to 
the degeneracy point, although B z ^ to avoid the environment-dominated regime. Fig.^Jshows examples of P(m, t). 
In (a) a measurement in the off-state, B x = 0, is shown, the other parameters are B z = 9.0 /LteV (corresponding to 
Sn g = 0.11) and \a\ 2 = 0.75. For comparison, Fig. 0Jb) shows a mesaurement where B x is on; the parameters are in 
this case B x = 1.3 fj,eV (which is the maximum value of B x assuming E™ ax = 0.2Eq), B z = 15 /LteV (corresponding 
to 8n g — 0.19) and \a\ 2 = 0.50. We see that after a while, P develops a two peak structure. This is interpreted as 
follows. At a given time, the total current is a superposition of two different currents flowing in the SET, with weights 
given by the amplitudes, a and f3, of the charge states in the qubit state. For increasing t, the peaks move towards 
higher m values since more electrons have tunneled through the SET. During this process, the peaks widen, their 
magnitude decreases and the distance (in m) between the peaks increases. The ratio between the two peak values are 
2.96 in Fig.Qfa) and 1.10 — 1.20 (the value decreases with time) in Fig.^fb) to be compared to the ratios |a| 2 /|/3| 2 = 3 
and 1 respectively. Fig. Ufa) corresponds to our proposed measurement situation - the off-state - and we see that 
here the peaks are well separated after a while and the ratio of their peak values stay close to |a| 2 /|/3| 2 for a very 
long time. Thus allowing for a slow accurate measurement of the charge. In Fig. 0Jb) on the other hand, the peaks 
are less well separated and the ratio agrees less well with |a| 2 /|/3| 2 . Eventually the valley between the peaks will fill 
out and one broad peak forms due to mixing between the charge states - a plateau indicating this process is clearly 
seen in (b). To perform a measurement of the charge, the peaks must have moved apart to become well separated 
and not yet started to merge due to mixing. This gives a window in time for measuring the charge - in the off-state 
in Fig.^fa) this window is very large - the lower bound is T meas ~ 500 ns. The upper bound is not seen, however, we 
expect it to be considerably larger. 



V. DISCUSSION 



We here estimate the decoherence times for the circular array. To be specific we assume Co = Cj — 100C g , 
Ec — 5E™ ax ~ 100/xeV and T ~ 40 mK = 3 /xeV and, initially N = 2 - higher N are commented on below. 

Calculations are performed by tuning H ctr i{t) in time. During this process the SET is off, and we are close to the 
degeneracy point. Assuming \Sn g \ < 0.2, then \B Z \ < 0.16^ — 16 ^eV - in addition we have \B X \ < 0.013i?c = 
1.3 /ieV. There is also a lower bound on AE, since the system is assumed to be in the Hamiltonian-dominated regime. 
The typical time per operation of the qubit is r op ~ h/ AE ~ 10 ps assuming AE ~ 10/ieV. 

For the ohmic noise, we use i? ~ 50 SI, which gives ot\ ~ 3 • 10~ 7 . Assuming AE < 16 /ieV one finds r v ~ 100 /is 
and T re i ax <; 100 /is. The Hamiltonian regime is realised when AE > aik B T » 10~ 12 eV. 

For l/f noise, t v is only known for 9 = and = ir/2. We use these two cases to estimate t v , assuming that 
this gives the correct order of magnitude for general 9. Nakamura et a/iii measured the factor af?t B = {Ef? f B ) 2 / E^ 
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for the SCB and obtained otfj^ B ~ 10~ 6 . If the l//-noise is caused by background charge fluctuations then it is 

reasonable to assume that for the CA — 2Naf9^ B , since the CA has 2N islands instead of the single island in 

the SCB. This gives E\^ = 2Trha uJook B T ~ 2NE^afj^ B and, for the parameters above, we find ~ 4 • lO -6 ^. 

From Nakamura et a/ii we also take u>i r ~ 310 Hz. For = we find T re i ax — > oo and r v ~ 1 ns if £? 2 < 16 /xeV 
and for = 7r/2, we find T re i ax ~ 10 ns and t v ~ 5 ns if B x < 1.3 /ieV. This shows that the Hamiltonian-dominated 
regime is realised if AE 3> 0.5 fieV. 

After the calculations, the SET is turned on and the measurement is started by tuning B x to 0. The relaxation is 
now very slow: T re i ax — ► oo when 9 — > for both ohmic and 1// noise. This gives ample time to measure \a\ 2 without 
using an ultra-fast detector, c.f. Fig.^ The dephasing is very rapid, thus the quantum state of the qubit is destroyed 
in a short time but it is still possible to measure the charge. This is good enough as a read-out mesurement but not 
as part of an error correction protocol. 

For higher N, the ohmic decoherence times will decrease somewhat, due to an increase in ot\. The same is true for 
1/ / noise, since E\/f ~ \ N. 

We conclude that it is the 1/ / noise that limits the operation of the CA as a qubit - it leads to the decoherence 
time Tjecoh ~ 1 ns. The decoherence due to ohmic noise is much slower. Hence, in practice, the ohmic noise seems 
to be of little concern. If the typical time for a quantum operation is r op ~ 10 ps, then the 1// noise restricts one to 
N op ~ 100 operations, assuming B x and B z are restricted to values that realise the Hamiltonian-dominated regime. 
This is a severe restriction and once again underscores that 1// noise is a serious problem for Josephson junction 
charge qubits. One possible solution to this might be the echo-technique demonstrated by Nakamura et alAL. 
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